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ABSTRACT 

Our current understanding of the physical conditions in the inner regions of proto- 
planetary discs is becoming increasingly challenged by the more detailed observational 
and theoretical explorations. Calculation of dust temperature is one of the key features 
we strive to understand and a necessary step in image and flux reconstruction. We ex- 
plore coexistence of small (0.1/im radius) and big (2/im radius) dust grains can coexist 
at distances from the star where small grains would not survive without big grains 
shielding them from the direct starlight. The study required a high resolution radiative 
transfer calculation capable of resolving large temperature gradients and disc surface 
curvatures caused by dust sublimation. The calculation was also capable of resolving 
temperature inversion effect in big grains, where the maximum dust temperature is 
at visual optical depth of ry ~ 1.5. We also show disc images and spectra, with dis- 
entangled contributions from small and big grains. Big grains dominate the near IR 
flux, mainly because of the bright hot inner disc rim. Small grains populate almost the 
entire inner disc interior, but appear in the disc surface at distances 2.2 times larger 
than the closest distance of big grains from the star. Nevertheless, small grains can 
contribute to the image surface brightness at smaller radii because they are visible be- 
low the optically thin surface defined by stellar heating. Our calculations demonstrate 
that the sublimation temperature does not provide a unique boundary condition for 
radiative transfer models of optically thick discs. The source of this problem is the 
temperature inversion effect, which allows survival of optically thin configurations of 
big grains closer to the star than the inner radius of optically thick disc. Future at- 
tempts to derive more realistic multigrain inner disc models will require the numerical 
resolution shown in our study, especially if dust dynamics is considered were grains 
can travel through zones of local temperature maxima. 

Key words: radiative transfer - circumstellar matter - protoplanetary discs - stars: 
pre-main-sequence - stars: imaging. 



1 INTRODUCTION 

A pressing issue raised by the ongoing proliferation of ex- 
trasolar planet discoveries is identifying conditions under 
which planets form and develop. A general scenario invokes 
a circumstellar (protoplanetary) disc of dust and gas around 
young pre-main-sequence stars, where planets are expected 
to form. The advancements in observational and theoret- 
ical techniques in recent decades enable d extensive explo- 
ratio n of the evolution of these discs ([Williams Ciezal 
l201lh . One of the key issues in protoplanetary disc explo- 
ration is understanding the evolution of dust properties. A 
major observational contributor to this field are invest iga- 
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tions of dust properties deduced form dust infrared (IR) ex- 
cess in the spectral energy distribution (SEP ) (e.g. various 
recent results from Spitzer surveys, such as ISargent et al.l 
2009; Watson 2009: Juh asz et al.P 2010: Mc Clure et all 2010 : 
Oliveira et al. 2010; M anoi et al.t 2011: Oliveira et allboill 
etc. and references therein). However, the exact details of 
the spatial disc structure are difficult to determine from 
the SEDs alone due to the intrinsic mathematical de- 
^eneracy o f models reproducing only flux measurements 
(IVinkovic et al. 2003 ). This makes direct imaging of proto- 
planetary discs an indispensable tool in exploration of their 
properties, but it also poses difficult challenges f or theo- 
retical m odels (e.^. Millan-Gabet et al 2007; Wat son et al.l 
120071 : iDullemond Monnierll2010l : IWilliams Ciezall201lh . 

The innermost disc regions, within a few AU around 
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the central star, have attracted a special attention in 
this field of research, since we expect formation of ter- 
restrial type of planets in this disc zone. It is also a 
zone of dust sublimation, high gas densities and com- 



plex dust dynamics (e.g. Cieslal l2009l: Vinkovi3 
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Armitage 2oTl|). At the same time, it is a very difficult ob- 
servational targe t due to its sm all scale and proximity to 
the central star (|Millan-Gabet et a l 2007). A strong impe- 
tus to the theoretical exploration of inner disc regions was 
given by development of semi- analytical models inspired 
by the observed properties o f the near IR part of SEDs 
(jPullemond k Monniei] l2010h . This is a part of the spec- 
trum where the inner disc emits thanks to its most resilient 
dust grains that survive up to temperatures of 1500-2000K. 

In recent years there is a tendency to improve upon the 
existing models. One motivation for that comes from ob- 
servations showing various features that cannot be accom- 
mo dated by the existing simplifi ed and time-steady mod- 
els (|Dullemond Monnierl I2OIOI ) . Another motivation for 
improvements spurs from theoretical investigations into the 
inner disc radiative transfer, grain and gas composition and 
disc geometry. For example, truncation of the inner disc 
by d ust sublimation results in a curvature of inner disc 
rim (llsella Nattall2005l : iTannirkulam. Harries Monnierl 
I2OO 7') : multigrain dust models res ult in a c omplicated ge- 
ometries and density distributions (jKama, Min Dominik 
I2OO9I I: big grains are gray in the near IR and experi- 
ence maximum temperature within the dust c loud and 
not at the very surface exposed to the starlight (|Vinkovicl 
I2OO6I : iKama. Min Dominil] I2OO9I ): dust dynamics is in- 
fluenced by rad iative pressure and magnetic fields in com- 
plicated ways (|Vinkovicl l2009l : iTurner Carballido Sanol 
I2OIOI I: dust and gas are thermally decoupled in optically 
thin regions, which makes gas hotter than initially expected 
(|Thi. Woitke KamDll201lh : gas itself can be a significant 
source of near IR flux and/or optically thin dust survives 
closer to the star than the optically thick i nner disc rim 
(|Tannirkulam et aLllioOSI : iBenistv et aDl2010l l. 

In this series of papers, we will investigate various theo- 
retical and numerical improvements in modelling of the inner 
disc structure. We start with a study of improved numerical 
reconstruction of dust thermal structure under a condition of 
multigrain dust composition. We use passively heated discs 
as a necessary initial step in identifying the most basic radia- 
tive transfer effects. Our focus is on coexistence of small and 
big dust grains, because these two dust properties yield big 
differences in dust sublimation behaviour. When used alone, 
small grains sublimate at considerably larger distances from 
the star than big grains. However, when mixed together they 
can coexist at distances dictated by big grains, but quan- 
titative details are still not known. Another problem arises 
from a complicated behaviour of big grains' temperature un- 
der direct exposure to stellar heating. As mentioned above, 
temperature initially increases with optical depth, reaches 
its maximum at some optical depth, and then decreases as 
we move deeper into the optically thick disc interior. Since 
dust sublimation is typically used as a unique boundary con- 
dition for inner disc structure, we will explore how solid is 
that assumption under high resolution numerical reconstruc- 
tion of multigrain dust temperature. In section [2] we describe 
some caveats in numerical treatment of the problem; in sec- 
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Figure 1. Cross section of olivine dust grains toorsch ner et akl 
1995) of 0.1/im and 2/im radius used in modelling. The dashed 
lines represent emission, the dotted lines scattering and the solid 
lines total extinction cross section. 



tion[3]we described obtained temperature structure and disc 
images and spectra; and in section 2] we discuss implications 
of our findings. 



2 RADIATIVE TRANSFER MODELLING 
2.1 Description of the model 

The main focus of this work is on the high resolution radia- 
tive transfer and analysis of dusty disc surface. Hence, we 
simplify the problen i of vertical disc structure and use the 
model developed bv IShakura Sunvae where the 

gas and dust are well mixed and the number density of dust 
particles of type a in the cylindrical coordinate system {g, z) 



na.{Q,z) ^ NocQ ^exp(-/io^— 



(1) 



with parameters p — 2^ m — 1.25 and ho — 1800. The 
coordinates (^, z) used in our calculations are scaled with 
the smallest distance of dust from the star, which we call the 
inner disc radius Rin and it is located in the disc midplane. 
This distance is not known in advance and it is a part of the 
output from radiative transfer calculation. The outer disc 
radius is fixed to Rout — lOORin- The number density of dust 
grains a at ^ = 1 and 2; = 0, which is the location of radius 
Rin , i s equal to iYo; = n^fl, 0) . We use olivine dust opacities 
from Dorschner et al.l (|l995l l and two coexisting grain types 
a of 0.1/im and 2/im radius. Their cross section is shown in 
Fig. [1] The ratio of their number densities at (1,0) is fixed 
to A^o.i : N2 = 10^ : 1. The dust sublimation temperature 
is Tsuh — 1500K. The spect r al sha pe of the stellar radiation 
is taken from the iKuruc i (119941 ) model atmosphere for a 
T* = 10 OOOK star. 

Instead of fixing the absolute value of dust densities, 
our radiative transfer calculation works with a fixed optical 
depth along an arbitrary radial line. We use visual optical 
depth — 10^ at 0.55/im toward the star in the disc mid- 
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plane. This optical depth is a line integral of dust visual 
cross sections {(7^^ ^ and cr^^^) and dust number densities 
(rio -L and 77-2), but we do not know in advance the integral 
limits because our two dust types sublimate at different dis- 
tances from the star. We expect 2/im grains to survive closer 
to the star than O.l/xm grains because they are better emit- 
ters in the near IR (see Fig. [T]), where Tsub has its spectral 
maximum. This means that 2fim grains are cooled more ef- 
ficiently and exist at Rin, while the sublimation front for 
0.1 /xm is not known in advance. 

We will see in the end that in our optically thick exam- 
ple smaller grains survive very close to Rin thanks to bigger 
grains shielding them from the direct stellar radiation. This 
simplifies the optical depth integral 

Rout/R^r^ 

V= J {(^o.i,v'^o.i(Q^^) + ^2,v'^2{Q^^))Rindg (2) 
and we can derive the dust density scale using equation [1] 



iVo.i 



+ Cr2,V^0-l/^V (1 - R^n/Rout) Rz. 



(3) 



This equation shows that density scale is directly related 
to the distance scale Rin. Using values for our parameters 
in equation [3] gives A^o.i = 3.76 x l(f / [Rin / RQ)m~^ . Since 
we do not know Rin in advance, we also do not know the 
absolute density scale. 

To solve the radiative transfer problem we only need 
to specify Tsuh and T*, the optical depth r^, the spectral 
shape of dust cross sections and the shape of stellar spec- 
trum. All other properties can be expressed in dimensionless 
terms. This is a consequence of intri nsic degeneracy of radia- 
tive transfer equations, as shown bv llvezic Elitzurl (|l997l l 
in their derivation of general scaling properties of the radia- 
tive transfer problem for radiatively heated dust. Luminosity 
L* is irrelevant and it enters as a parameter only after we 
obtain Rin in units of stellar radii i?* from the output of 
radiative transfer calculation. Then we can scale our model 
according to the stellar radius derived from the luminosity: 
R. = RQ{TQlT.Y{U/LQf-\ 



2.2 Boundary condition at the inner rim 

The parameter that has been traditionally used as a bound- 
ary condition for the inner rim of dusty discs is Tsuh- The 
physical motivation for this was the expectation that the 
dust temperature of externally heated dusty clouds always 
monotonically decreases with the increasing distance from 
the heat source. This means that the maximum dust tem- 
perature is achieved at the closest distance to the star that 
a dust particle can reach before sublimating away. 

However, this concept was found incorrect in the case of 
externally heated optically thick clouds made of big (> l/xm) 
dust grains. Such grains have an ability to efficiently absorb 
the diffuse IR radiation originating from the clouds interior, 
which results in an increase of temperature within the dust 
cloud at the visual optical depth of ^ 1, instead of on the 
very surface exposed to the stellar heating. This process of 
temperature inversion is similar to the greenhouse effect and 
it was first discovered numerically in single grain size models 
(lDullemondl200i : llsella Nattall2005l : IVinkovic et alll2006l l 
and then proved analytically for multigrain dust mixtures 




Figure 2. Computational grid used in the final model. The panels 
show two spatial scales of the inner disc region in cylindrical coor- 
dinates z). Thick dashed lines are radial visual optical depths 
of = 0.1 and ry = 1. A section of the grid with increased 
resolution traces the surface of disc populated with O.l/xm grains, 
while 2)U,m grains exist at all grid vertices. 



(jVinkovid [iooel l. The analytical analysis showed that big 
grains dominate the cloud surface, while smaller grains can 
exist immediately behind the zone of temperature inversion. 
The analysis also showed that the effect exists for both types 
of radiative transfer boundary conditions: (1) a constant ex- 
ternal flux heating the cloud, with no limits on the dust 
temperature, and (2) a fixed maximum dust temperature, 
corresponding to dust sublimation. 

The inner regions of protoplanetary discs provide ideal 
conditions for emergence of the temperature inversion phe- 
nomenon. Accretion, which is ubiquitous to protoplanetary 
discs, constantly resupplies the inner regions with big grains 
that grow in the disc. Hence, the optically thick discs should 
be in a permanent process of minimization of their inner disc 
radius, with big grains populating the inner disc surface and 
dictating the radiative transfer. A complication arises when 
we try to deal with the spatial scale of the disc surface be- 
cause the surface is defined by the optical depth of < 1 as 
a zone where the most of stellar radiation is absorbed. This 
optically thin zone is, therefore, not necessarily geometri- 
cally thin. Instead, it could extend much closer to the star, 
creating a large optically thin dusty zone spreading over 
radii smaller then the inner radius of optically thick disc. 
Since dust can exist only up to its sublimation temperature, 
IVinkovi(3 (|2006l l proposed using two inner radii that discrim- 
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inate between these two coexisting, optically thin and thick, 
disc zones: 



Tsuh 



(4) 



where ^ = 2 is used for the optically thick 
(|Dullemond. Dominik Nattal l200l | l and ^ ^ 1.2 for the 
optically thin inner radius (|Vinkovic et al.ll2006l ). Note that 
here we assume that the disc opacity is dominated by big 
grains that have gray opacity in the near IR. Smaller grains 
are inefficient emitters of the near IR radiation and, there- 
fore, overheat and sublimate away at these distances form 
the star. Smaller grains can survive either at larger distances 
fro m the star or hidden in the disc interior behind gray dust. 

iKama, Min Domin ik (2009) showed these effects nu- 
merically, using discs with various dust compositions and 
surface densities. The optically thick radius in their models 
dominated by big grains are in agreemenlQ with equation 
m The optically thin zones in their models are dictated by 
predefined power law density profiles and, therefore, do not 
extend to the t heoretical minimum distance t o the star. 

Models bv lKama, Min &: Dominic (|2009l l demonstrate 
the boundary condition problem that we face when the spa- 
tial dust density distribution of the optically thin zone is 
free to vary. In order to simplify the problem and make it 
manageable, in this paper we do not explore properties of 
the optically thin zone. Our goal is limited to the explo- 
ration of temperature inversion structure in 2D and survival 
of small grains in the disc. We achieve this by fixing the 
maximum dust temperature anywhere in the disc to Tsuh- 
This means that we do not know the dust temperature Tin 
at = 1, z = 0), but we expect it to be smaller than Tsuh- 

Such a boundary condition is easy to postulate, but it 
creates substantial numerical challenges, especially in a case 
of high resolution temperature calculation like ours. Our nu- 
merical approach was semi-autonomous, where the compu- 
tational grid updates and Tin corrections where checked by 
hand for precision and error control. We know in advance 
that 2/im grains in our model should survive at all grid 
points. On the other hand, O.l/xm grains must be removed 
from the surface of inner disc where these grains overheat. 
They survive only in the interior shielded by 2/im grains, 
but their exact location of survival is not known in advance. 

The solution is looked for by iterating between radia- 
tive transfer calculations and dust removal. These iterations 
were combined with corrections by hand of Tin- In addition 
to that, we updated the computational grid structure after 
each dust removal because our grid is irregular and auto- 
matically traces spatial and optical depth gradients. Even 
though the majority of this iterative process is automatic, 
we had to check each step by hand to avoid errors and grid 
imperfections. Nonetheless, small deviations from Tsuh are 
very difficult to suppress, but we managed to keep the error 
within 3% (45K for Tsuh = 1500i^^). 



^ Small deviations from equation 3] are visible in their models 
because they use dependence of dust sublimation on gas density, 
which slightly increases dust sublimation temperature at the inner 
disc rim. 



2.3 Radiative transfer calculation 

Numerical radiative transfer calculations were performed 
with the code LELUYA (www.leluya.org). It solves the in- 
tegral equation of the formal solution of radiative transfer 
with axially symmetric dust configurations, including dust 
scattering, absorption, and thermal emission. The solution 
is based on a long-characteristics approach to the direct 
meth od of solving the matrix version of the integral equa- 
tion teurucd ll969'). The equations are solved on a highly 
unstructured triangular self-adaptive grid that traces simul- 
taneously both the density gradients and the optical depth 
gradients over many orders of magnitude in spatial and op- 
tical depth space. The code is parallelized and written in 
C. Th eoretical and computational details are explained by 
IVinko vic (2003). 

The key for a high resolution temperature structure is 
the grid used in our calculations. The code creates a grid 
by integrating optical depth along lines that cut through 
the largest optical depth gradients. In our case these are 
radial lines and vertical lines parallel to the symmetry axis. 
The grid vertices are chosen on these lines according to the 
following equation: 



At AL 



(5) 



where Ar is the optical depth step, AL is the spatial step, 
where both are calculated relative to the last created grid 
vertex, and r and L are the optical depth and the spatial 
distance from the starting point of integration. If used care- 
fully, with starting points in low density areas and at dust 
sublimation surfaces, this equation enforces higher grid res- 
olution in zones of optical depth gradients. Vertices in op- 
tically deep interior are distributed more sparsely because 
the radiation field is slowly varying. In addition to this au- 
tomatic grid generation, we also inspected the grid by eye 
and added extra grid vertices where we considered necessary. 
The final grid is shown in Fig. [2l 

Each grid vertex contains an angular grid defining direc- 
tions of integration of radiative transfer equation. It would 
be simple and fast to use always the same type of angular 
grid, but vertices on the disc surface are surrounded by an 
anisotropic radiation field that changes its geometry from 
vertex to vertex. Therefore, we use a method where a unit 
sphere is split self-adaptively into small spherical triangles 
according to the number of vertices visible through a trian- 
gle in 3D space. This enables increased resolution in direc- 
tions of radiation field gradients. The surfaces of spherical 
triangles are used as corresponding statistical weights in the 
discretized angular integration. 



3 RESULTS 

We obtained the inner disc radius of Rin = 44.72i?*, 
which agrees with equation [4l for the case of optically thick 
dusty disc edge. In real units it corresponds to Rin — 
O.O687(L*/L0)°-^AU for a 10,000K star. Typical luminosi- 
ties of Herbig Ae stars are L* = 20 — lOOL©, hence, typ- 
ical values for the inner optically thick disc radius are 
R,^ = 0.31 - 0.69AU. 
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Figure 3. Temperature of dust grains shown at two spatial scales in cylindrical coordinates for two coexisting grain radii (indicated in 
images). The dashed line indicates the radial optical depth of ry = 0.1, while the solid line is ry = 1. Note how 0.1/im grains survive in 
the optically thin disc surface only at radial distances of ^ > 2.2, but in the disc interior they exist very close to Rin- The dotted lines 
in the most upper panel indicate radial lines of 88° and 89° inclination angle used in Fig. |4] and [5] The dotted line in the third from the 
top panel shows the vertical optical depth of ry = 1. 
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Figure 4. Temperature profiles along radial visual optical depth 
lines for different inclination angles relative to the symmetry 
axis (z-axis): 90° (midplane), 89° and 88° (shown in Fig. [3] as 
dotted lines in the top panel). The solid line indicates the tem- 
perature of 2/im grains and the dashed line of O.l/xm grains. The 
error bar marks 3% error tolerance on the dust sublimation tem- 
perature of 1500K. 
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Figure 5. Same as Fig. (4] except that the temperature is shown 
relative to the spatial scale R = {q^ + z^)^-^ along radial lines. 
Steep decline of 0.1/i,m grain temperature explains increased grid 
resolution in the sublimation zone of 0.1/im grains (see Fig. O. 



3.1 Temperature structure 

The dust temperature distribution is shown in Fig. [3] for 
both 2/im and O.l/xm grains. The optically thin disc surface 
is populated only by big grains at distances ^ < 2.2 from 
the star. However, smaller grains survive in the optically 
thick disc interior, where they are heated by diffuse infrared 
radiation and protected from the direct stellar radiation. 
This enables smaller grains to survive almost as close as 
Rin from the star. 

The temperature inversion effect is reproduced by our 
model, with the temperature maximum at ^ 1.5 in the 
optical depth space and at ^ 1.005i?in from the star. As a 
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Figure 6. Map of temperature along the line of sight for differ- 
ent inclination angles. The images show temperatures at optical 
depth r = 1 for 2.2/im wavelength. The left side of images shows 
temperatures of 0.1/im dust grains and the right side shows tem- 
peratures of 2/im dust grains. 
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Figure 7. Spectral energy distribution of the model for three 
different disc inclination angles: 0° (face on), 45° and 70°. 




0.1 1 10 100 

I (iLim) 



Figure 8. Wavelength variation of the relative contribution of 
stellar, emission and scattering flux to the total flux for 45° disc 
inclination. Other inclinations have very similar relative contri- 
butions. Note the dominance of scattering in the near IR due to 
big 2fim grains (see also dust cross sections in Fig. [T}. 

result of this effect, the temperature decreases very slowly 
with the optical depth and the zone of high dust temperature 
is not confined to a geometrically thin disc edge. With high 
temperature we mea n temperatures h igh enough to quickly 
anneal the silicates (jmil et al.ll200lh . which is >950K. In 
Fig. [3] we can see that the disc temperature is above 950K 
between Rin and ^ l.lRin, which is geometrically a signifi- 
cant part of the inner disc. 

A more detailed insight into the temperature structure 
is given in Fig. [Hand Fig.O They show cuts through the disc 
along radial lines of different angles 6 relative to the symme- 
try axis (z-axis). Fig. [H displays the dust temperature rela- 
tive to the radial optical depth. In the midplane (0 — 90°) 
temperature stays above the annealing limit of 950K almost 
until TV ^ 10^. Small grains sublimate at rv^ ~ 4 — 6 and 
their temperature drops faster then the temperature of big- 
ger grains, such that both grain sizes reach the same temper- 
ature at r > 10. This is i n agreement wit h the qualitative 
theoretical predictions by IVinkovid (|2006l ). The cut alone: 
= 88° shows a slower decline of 0.1/im grain temperature. 



Figure 9. Wavelength variation of the relative contribution of 
0.1/im and 2/im grains to the total observed diffuse flux for 45° 
disc inclination. Other inclinations have very similar relative con- 
tributions. Note how 2/im grains dominated the near IR because 
the opacity of 0.1/im grains drops fast in this wavelength region 
(see Fig. [T}. 



but this is because this cut goes through the disc region 
where the largest temperature gradient is perpendicular to 
the disc plane (see the upper dotted line in the top panel in 
Fig.E]). 

Fig. [5] displays the same radial cuts, but relative to the 
spatial distance from the star. Here we can see that the dust 
temperature in the midplane stays above 950K almost up to 
l.lRin- Small grains sublimate at distances from 1.008i?in 
(midplane) to 1.045i?in (for = 88°). The temperature gra- 
dient for O.l/xm grains is very steep between the point of 
Tsuh and the point of thermalization with 2/im grains. This 
explains why the computational grid has such a high resolu- 
tion in the region of O.l/xm dust sublimation (see Fig. [2]). The 
spatial step in this region goes as low as AL = 3.7x 10 ~^i?in. 

The temperature of 2/im grains is slightly increased at 
the location of sublimation of O.l/xm grains. The effect is 
small because the near IR emissivity of 0.1/im grains is too 
small to provide a significant contribution to the local dif- 
fuse radiation. This would not be the case if the number 
density of bigger grains were reduced to the levels where the 
diffuse radiation is dominated by smaller grains. However, 
this would also mean that bigger grains do not dominate 
the optical depth and, in turn, cannot provide shielding to 
smaller grains. This would result in a larger inner disc radius 
dictated by the opacity of smaller grains. 

Another way of observing the disc temperature struc- 
ture is to focus on how deep we can see into the disc. For 
example, the dotted line in the third from the top panel 
in Fig. [3] shows the optical depth of = 1 when the disc 
is viewed face on. This is approximately the depth of ob- 
servable disc at 0.55/xm wavelength. Fig. [6] shows the dust 
temperature distribution at the optical depth of r = 1 at 
2.2/im wavelength for different inclination angles. It reveals 
that we can see deep enough into the disc to sample both 



8 D. Vinkovic 




-4 -2 2 4 

X (RJ 



4 -2 2 4 

X (RJ 



Figure 10. Model images at 2.2/im (left column) and 10/im (right column) wavelength for different inclination angles. Images at 
wavelengths < Sfim show brightness behaviour similar to this example. A faint brightness ring in images at ^ 2Rin distance from the 
centre is a feature caused by 0.1/im grains in the disc surface. The flux through each image pixel is scaled with the bolometric flux -F^oZ- 
Contributions from 0.1/xm and 2fim dust grains to the model images is shown in Fig. 1111 
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Figure 11. Contribution from O.l/xm and 2/im dust grains to the model images at 2.2/im (left column) and 10/im (right column) 
wavelength for different inclination angles. The left side of images shows contribution from 0.1 fim dust grains and the right side from 
2fim dust grains. The flux through each image pixel is scaled with the bolometric flux -F^oZ- For the complete images merging both 
contributions see Fig. 1101 
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Figure 12. Model image at 2.2/im (left column) and 10/im (right column) wavelength and 45° inclination separated into its dust scattering 
and emission component. The left side of images shows contribution from 0.1/im dust grains and the right side shows contribution from 
2/im dust grains. The total image is shown in Fig. 1101 



dust types even in the regions where smaller grains do not 
populate the disc surface. 

3.2 Spectra and images 

The SED of protoplanetary discs traces the dust emission 
and scattering at different temperatures and optical depths. 
The inner disc region experiences the highest temperatures 
and it should be detectable at near IR and mid IR wave- 
lengths. Large optical depths prevent photons of these wave- 
lengths from the deep disc interior to be detected; hence, the 
optically thin disc surface dominates the flux. Fig. [71 shows 
SEDs of the model viewed face on (inclination angle 6 — 0°) 
and inclined {0 = 45° and 6* = 70°). The SEDs do not differ 
at short wavelengths where the stellar flux dominates. They 
also do not differ at the far IR wavelengths where the disc 
becomes optically thin to its own emission. 

The relative contribution of stellar flux, dust emission 
and scattering to the total SED is shown in Fig. [8] for the 
face on disc model. Inclination decreases scattering contri- 



bution and increases the contribution of dust emission, but 
the overall wavelength variation is similar to this face on 
example. The competition between dust emission and scat- 
tering for the dominance in the infrared part of spectrum 
is closely matched with their interplay in the dust opacity 
shown in Fig. [T] The dominance of big grains in the disc 
surface leads to strong contribution of dust scattering in the 
near and mid IR. The details depend on the chosen dust 
properties, but the overall trend indicates that scattering is 
a significant contributor to the SED. 

The near and mid IR part of the SED decrease approxi- 
mately as cos 6 for small or intermediate angles. The surface 
is geometrically thin and cos^ comes from the geometrical 
projection of the emitting surface. For larger inclination 
angles the SED decreases slower than cos^, since the line 
of view is grazing the disc surface and its geometrical 
thickness is not negligible any more. A somewhat different 
SED behaviour is visible at 2 — 3/im where the flux is 
almost independent of inclination for small angles and 
then decreases only slightly with larger angles. These are 
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wavelengths where the inner disc edge emits at Tgub tem- 
perature. The curvature and vertical height of the inner disc 
edge enable its emission to be almost inse nsitive to 0, un- 
til t he incl ination becomes large (Isella. Testi Nattal 



20061: 



.Tannirkulam, Harries Monnierl l2007l : 

Thi. Woitke KamDll20m 

Model images shown in Fig. [10] demonstrate this. The 
inner disc edge is the brightest feature in 2.2/im images and 
inclination exposes its curvature. At lO/xm the disc edge is 
not the main feature and the rest of the disc surface dic- 
tates the SED. Disentangled contribution from both grains 
sizes to the disc images is shown in Fig. 1111 The bright- 
ness of the inner disc edge is dominated by big grains, but 
smaller grains also contribute despite being positioned op- 
tically deeper within the disc. A faint ring visible at ^ 2.2 
distance from the image centre is the location where smaller 
grains show up in the disc surface. 

The near IR images and SED are obviously highly sen- 
sitive to the curvature of the inner disc edge. This is why its 
height and shape has been of great interest to astrono mers 
(|Millan-Gabet et alllioOTl i lDullemond &: Monnierll20ld ). But 
the rest of the disc surface also contributes to the image 
flux. In Fig. [12] we show how dust scattering and emission 
differ in their contribution to the disc surface brightness at 
2.2/im and 10/im wavelength. The thermal emission is mod- 
ulated by the dust emissivity and reddened by dust layers 
along the line of sight. Fig. [6] shows the temperature distri- 
bution across the 2.2/im image at r = 1 optical depth, which 
approximately shows how deep we see into the disc. From 
Fig. [12] we see that temperatures of 2/im grains contribut- 
ing the most to the thermal emission at near IR wavelengths 
are limited to small radii close to the inner disc edge, while 
0.1/im grains contribute over larger area, but with smaller 
surface brightness due to their smaller emissivity. In the mid 
IR both grains have thermal emission spread over a larger 
area. Since O.l/xm grains appear in the disc surface layer at 
Q ~ 2.2, we can see a brightness increase in images at this 
distance from the centre. However, this increase is small be- 
cause the emissivity of O.l/xm grains is small compared to 
2/im grains. This is why big grains dominate the SED and 
images in the near IR even though smaller grains have high 
temperatures over a larger portion of the disc surface. 

While the emission is confined to small radii where the 
dust is the hottest, scatter ing follows s urface dust density 
and extends further away. IPinte et "al] (l2008l l showed how 
scattering is an important component in the interpretation 
of the near IR imaging data. We need to emphasize here that 
our images are calculated with isotropic scattering function, 
which is a rough approximation for big grains. A more realis- 
tic scattering function would yield larger image asymmetries 
with inclination. 



4 DISCUSSION AND CONCLUSIONS 

In this paper, we have explored the temperature structure 
of the hottest dusty regions of protoplanetary discs. In this 
the most inner part of the disc, the interplay of dust density, 
optical depth and multigrain dust composition creates an 
intricate disc surface structure shaped by nonuniform subli- 
mation of individual dust species. Our study was limited to 
two distinctively different dust grain radii of 0.1/im (small) 



and 2/im (big). These grains exhibit very different sublima- 
tion behaviour when used individually in single grain mod- 
els, with big grains surviving much closer to the star than 
small grains. However, we demonstrated a significant de- 
crease of sublimation distance of small grains when they are 
mixed with big grains under conditions of big grains domi- 
nating the opacity. Small grains survive behind big grains' 
visual optical depth of tv 4 — 6, which translates into 
only a few percentages larger sublimation distance than big 
grains. We argue that the disc accretion constantly resup- 
plies the inner disc with big grains, which enables them to 
dominate the opacity and sustain conditions described in 
our paper. There are also observational confirmations of this 
inner disc property; spectra indicate larg e fractions of big 
grain s in the surface of inner discs (e.g. Ivan Boekel et al 



grain s m tne suriace oi inner discs (e.g. [ van Jjoekel et al 
2004) . I2OO5I : [Sargent et al.l l2009l : IWatsonll2009l : ljuhasz e~ 



while the observed sizes of inner discs are consis- 
tent with big grains dic tating the inner disc radius Rin 
(jMillan-Gabet et alll2007 h. 

Our study required very high self-adaptive computa- 
tional grid resolution (almost 10~^ Rin) in order to cor- 
rectly trace dust sublimation and temperature gradients. 
Computations have to resolve two important tempera- 
ture gradients. One is the temperature inversion experi- 
enced by big grains, where the temperature reaches its 
maximum ai rv ^ 1.5 instead of the very surface ex- 
posed to the stellar radiation. The other comes from the 
tendency of small grains to quickly thermalize with big 
grains as the optical depth is increasing, which results 
in a steep temperature decline for small grains. Prior 
similar studies wer e not resolving these effects in multi- 
grain dust models (iTannirkulam, Harries Monnierl l2007l : 



iKama, Min DominikI 2009h , which signifies the impor- 
tance of high resolution radiative transfer calculations. 

Implications of the temperature inversion effect, where 
the dust temperature at Rin is lower than sublima tion, have 
received modest attention in the lit e rature so far (IVinkovic 
| 2006l : iKama. Min Dominil3 l2009l : iDullemond Monnier 
2010; . The main consequence of this effect is that big grains 
can survive closer to the star than optically thick part of 
the disc. For theorists this implies that currently used dusty 
disc models are incomplete, while developing a more realis- 
tic model will require scrutinizing the dust dynamics coupled 
with radiative transfer. For observers this implies that op- 
tically thin dust might be detectable within the currently 
used optically thick inner disc radius. Recent detection of 
disc brightness within Rin could be the first indication that 
such a dust structure exists, bu t it has to be disentan- 
gled from possible gas emission ([Tannirkulam et al ] I2OO8I : 
iBenistv et"aDl2010h . 

We also resorted to a simplified boundary condition 
where the dust density drops to zero at Rin, simply be- 
cause we currently do not have any theoretical handle on 
the density structure of optically thin dust within Rin- Its 
structure can be highly infiuenced by dust dynamics due to 
other forces than gravity and gas drag. Radiation pressure 
force would act upon all dust grains in this optically thin 
disc zone. Magnetic fields at these distances could be strong 
enough to interact with charged dust particles. Charging is 
expected because of high gas temperatures and exposure to 
UV and X-ray radiation ([Turner, Carballido San 3 I2OIOI : 
IPedersen &: Ines Gomez de Castroll201lh . This leads to cou- 
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pling of grain settli ng to the distribution of the magnetoro- 
tational turbulence ([Turner, Carbalhdo &: .Sario|[2010).High 
gas densities also provide conditions for photophoresis on 
big grains, which moves particles outward (Moudens 2011). 
Moreover, these dynamical processes could produce optical 
depth variability within Rin, which would immediately in- 
fluence the optically thick region where the dust at the sub- 
limation zone can overheat and sublimate. Since this subli- 
mation zone is located within a disc, we do not know how 
optical depth changes caused by dust sublimation would in- 
fluence the overall disc structure. All these issues make fur- 
ther work on the tempera ture structure of inner disc region 
a very challenging task ( Dullemond MonnierlbolO ) that 
will require computation of individual trajectories of dust 
grains coupled with radiative transfer. Such models are still 
under development for studying the outer, dynamically and 
thermally simpler disc regions (e.g. Charnoz et al. 2011). 

We have also analyzed how differences in the distribu- 
tion of small and big grains influence the disc spectrum and 
images. Big grains dominate the near IR spectrum and im- 
age flux, since small grains cannot survive in the disc sur- 
face up to distances > 2.2Rin- Nonetheless, small grains 
contribute signiflcantly to the image surface brightness, es- 
pecially at scales several times larger than Rin. Hence, in- 
voking grain growth as an interpretation for a larger fraction 
of big grains in the inner part of protoplanetary discs has to 
take into account the dominance of big grains in the near IR 
opacity of dust mixture and the spatial separation of dust 
grains by dust sublimation. Notice that in our model we 
used the same ratio of big to small grains through the entire 
disc, but, nevertheless, their fractional contribution to the 
SED and images is highly wavelength dependent. 

Dust separation is also caused by radiation pressure 
force, which we will explore in the next paper. Small grains 
are more easily eroded from the surface of optically thick 
discs than bigge r ^ains due to the rad iation pres sure from 
stellar radiation (jXakeuchi Linl l2003V Moreover, IVinkov"i3 
(I2OO9I I showed how big grains experience additional radia- 
tion pressure force due to the near IR radiation from the 
hot disc dust, which pulls bigger grains out of the disc until 
the stellar radiation pressure force overcomes the gas drag 
force and blows grains outward. This process is very fast 
(within orbital time) unless the local gas density is high 
enough to dominate over radiation pressure force. Another 
important effect is corrections to the vertical size of the in- 
ner disc. The inner disc edge is hotter than the standard 
disc structure, which results in vertic a l expansion of the disc 
(iDullemond. Dominik k Nattal [2OO1I : iThi. Woitke k Kamd 
l201lh . As long as the disc does not expand so much to cause 
self-shadowing, we expect its temperature structure to be 
qualitatively similar to our model, except for geometrical 
scaling of the vertical disc size. Another important thermal 
effect comes from viscous heating, which increases the tem- 
perature in the midplane. However, here we again expect 
the thermal structure above the midplane to be dominated 
by effects described in o ur model, since t his region is domi- 
nated by stellar heating (jD'Alessio et al.ll2006l ) . Overall, the 
exact dynamical stability of dusty disc structure taking into 
account all these effects, as well as optically thin dust within 
Rin^ is still unknown and it requires additional research. 
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